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Abstract 

It is shown that due to memory effects the complex behaviour of components in a 
stochastic system can be transmitted to macroscopic evolution of the system as a whole. 
Within the Markov approximation widely using in ordinary statistical mechanics, memory 
effects are neglected. As a result, a time-scale separation between the macroscopic and 
the microscopic level of description exists, the macroscopic differential picture is no a 
consequence of microscopic non-differentiable dynamics. On the other hand, the presence 
of complete memory in a system means that all its components have the same behaviour. 
If the memory function has no characteristic time scales, the correct description of the 
macroscopic evolution of such systems have to be in terms of the fractional calculus. 

PACS number(s): 05.40.-a, 05.45.-a, 05.60.-k 



1 Introduction 



The connection between microscopic dynamics of components in stochastic sys- 
tems and macroscopic description of their behaviour as a whole is very attrac- 
tive in statistical physics [|I], @, 0]. From the point of view of a single trajec- 
tory, the path of a Brownian particle is continuous everywhere, but is nowhere 
differentiate ||]. Such a path is not described by an analytical function. In 
the theory of Brownian motion formulated by Langevin[[5] the velocity of the 
Brownian particle was proved to be discontinuous ||. The differentiate nature 
of the macroscopic picture of Brownian motion is due to the key role of the 
Central Limit Theorem (fluctuations of the microscopic quantities are indepen- 
dent of each other). This means that the microscopic and macroscopic levels 
of description of the process are separated in the time scale, and memory of 
the non-differentiable character peculiar to the microscopic dynamics is lost in 
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the long time limit. Consequently, the results of observing the motion of an 
ensemble of trajectories can be predicted by means of theoretical prescriptions 
based on ordinary mathematical procedures proceeding from the differentiabil- 
ity assumption. When the condition of time-scale separation is not available, 
the non-differentiable nature of the microscopic dynamics can be transmitted to 
the macroscopic level fl7j. In the present paper we show that the key to this un- 
derstanding is memory effects in stochastic systems. It is due to memory effects 
that the macroscopic behaviour of stochatic systems contains a manifestation 
of microscopic dynamics. 

The outline of the paper is as follows. Section g is devoted to the behaviour 
of relaxation in the physical systems without and with complete memory. In 
Section [3] we use the generalized (in the terms of the memory effects) Langevin 
and Kramers-Moyal equations to shown that, by means of the memory func- 
tion, the non-differentiable nature of microscopic dynamics of system compo- 
nents can be transmitted to the macroscopic level of description in the form 
of fractional derivatives. Note that the memory effects can induce orderedness 
of macroscopic processes in stochatic systems. In Section f| we briefly con- 
sider the criterion of the relative degree of order in the systems. Section |5] is 
devoted to the Time Fractional Diffusion- Wave equation concerning the case 
when chaos and order coexist. Finally we compare the process described by the 
latter equation with the fractional Brownian motion. 

2 Classification of Memory Effects 

We start from some classification of memory effects. It is based on mathematical 
properties of the corresponding memory function. Let us consider the integro- 
differential equation 



where / is the quantity of interest, K the memory kernel, and A the parameter. 
The equation is a typical non-Markovian equation obtained in studying the 
physical systems coupled to an environment, with environmental degrees of 
freedom being averaged. The parameter A can be regarded as the strength of 
the perturbation induced by the environment of the system. 

For a system without memory (ideal Markov system), the time dependence 
of the memory function K(t — t') is of the form 




(i) 



K(t - t') = S(t - t'), 



(2) 
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where 5(t — t') is the Dirac (5-function. The absence of memory means that the 
convolution function, J(t) = Jq K{t — t') f(t') dt', is defined by f(t) at the only 
instant t. Substituting ([2]) into (|TJ) we obtain 

Equation (|3|) has an exponential solution. If memory effects are introduced 
into the system the (5-function in (0) turns into a bell-shaped function, with the 
width determining an interval r during which f{t) has an effect on the function 
J. 

In the literature, wide use is made of the Markov approximation, which 
replaces the equation ([I]) with the rate equation 

<m = -^f mdt >) m . (4) 

The Van Hove limit [|8f implies that the limit A — > 0, t — > oo is approached to 
in such a way that the product \H is kept constant. Using the Van Hove limit 
makes it possible to replace the time convolution in (|T|) with 

f> - -AVW, (5) 

where 

r = fdtK(t). 

The limit A — > implies that the coupling of the system to the environment is 
weak, while the limit t — > oo means that the observation time is much larger 
than the characteristic time scale r. 

On the other hand, in the systems having ideal complete memory the func- 
tion J is formed over all the course of the action of the quantity f(t') up to the 
instant t with the weight K{t — t') = {1,0 < t' < t;0,t' > t] (step function). 
In this case the equation ([]]) is transformed into 

f> = -A7>^. (6) 

It has the unique solution, f{t) = /(O)cos(At), which does not decrease at 
t — > oo in contrast to (0). 

Relation ([I]) written in the time domain is not always convenient because of 
the convolution (integral over t'). This can be eliminated by using the Laplace 
transformation 

f(s) = C[f(t)} = J f(t) e- st dt, f(t) = —J f(s) e st ds. 

JU /7T7, J— too 
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(7) 



where the initial condition is taken into account. The Laplace transform of the 
kernel (0), which corresponds to the absence of memory, yields the constant, 
K(s) = 1. For ideal memory, we obtain K(s) = 1/s. Thus, as the ideal 
complete memory appears in the system the constant kernel is replaced by the 
hyperbolic one. It is logical to infer that the fractional integration of the order 
z/, < v < 1, will interpolate the memory function between the 5- function and 
step function. Systems with such a memory function occupy an intermediate 
position between the two limiting cases and are described in [|9f. They have 
complete but not ideal memory. This means that the memory manifests itself 
within all the interval (0, t) but not at each instant t' . Such a memory function 
has no characteristic time scale, i. e. the Markov approximation is inapplicable. 

Assume that memory holds only at the points of a Cantor set. The problem 
lies in finding the Laplace transform of a step memory function. To construct a 
Cantor set, we first choose the entire time interval of the length T and remove 
the central part of the interval leaving two intervals of the length £T (where 
£ < 1/2). Obviously, to avoid the loss of the integral memory, the heights 
of the two resulting bars must be increased to the value (2£T) _1 . In the next 
stage, each remaining interval of the length £T is subjected to the same division 
process. In each subsequent stage n, this contraction procedure is performed 
for the 2 n_1 stages obtained in the preceding stage. One can easily see that the 
memory function K n (t) is represented by a set of 2 n bars of the height 1/(2£)"T 
and of the width £ n T. The Laplace image of K n (t) is written as 



For n ^> 1 we have |sT£ n | <C 1. As the special investigations show 
(Chapter 5), the limiting value of the function K n (s), when the number 
of divisions generating the Cantor set tends to infinity, becomes 



where v = ln2/ln(l/£) is the fractal dimension and q(\n(sT)) is the periodical 
function with the period ln£. From the physical point of view the fractal 
dimension v informs us about the ralative amount of states being conserved 
in the process of interaction and represents a quantitative measure of memory 
effects. It is clear from the analysis of the limiting case. For an empty Cantor set 
[y = 0) the dependence (0) reduces to the constant, corresponding to the entire 



K n (s) 



1 - exp(-sT£ n ) n-i 1 + exp(-z£ k ) 
11 2 



z = (l-()Ts. (8) 



K(s) = (sT)- u q(\n(sT)) 



(9) 
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absence of memory (in the ^-domain it consists of two delta-functions located 
on the edges of the Cantor set). The limiting value of the similarity parameter 
£ = 1/2 yields the dimension v = 1, which corresponds to the complete memory. 
The result (|9j) is correct for any sT from the interval 

i/(i-0 < \ sT \ < 00 

or for the time variable located in the interval 

0<t/T<l. 

Averaging the function q(ln(sT)) in (|9|) over the period ln£ (see details in 
H) and taking the inverse Laplace transform of J(s), we obtain the temporal 
fractional integral 



where r(z/) is the Gamma function (see Appendix). The averaging procedure 
converts the discrete fractal density p n = l/(2£) n T to its continuous value, 
oc Now the kernel K(t - t') in J|) becomes 

- *') ~ j^y(* - 

indicating the presence of the fractional integral. The integral representation 
of equation (fj) is equivalent to a differential equation of the fractional order. 
The above result (|9|) can be generalized to the case of an arbitrary number 
j of elementary blocks participating in the construction of a Cantor set. It 
turns out|9j that the result (j9[) remains true given that v = lnj/ln(l/£). The 
transition from the regular Cantor sets to the case, when the parameter £ is 
random in each stage of constructing a Cantor set, leads to the same result. 

Thus, the systems with residual memory have K(s) — s _i/ , where the expo- 
nent value, < v < 1, determines the extent of memory preservation. Substi- 
tuting it into (0) (let all constants be contained in A 2 ) and using the inverse 
Laplace transformation, the solution of equation (JXJ) takes the form: 

f(t) = /(0) E V+1 (-\H V+1 ) , 

where E a (z) = EkLoZ k /T(ak + 1) is the one-parameter Mittag-LefHer 
function |TDJ. In the particular cases v = and v = 1 we have 

Ei(-z) = e~ z , E 2 (-z 2 ) = cosz. 

For < v < 1 the quantity f(t) has an algebraic decay as t — ► oo. There- 
fore memory effects can essentially change the character of relaxation. The 
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above constructions are not only pure theoretical, but reflect the experimental 
situation. In this connection it should be mentioned that, for example, the 
relaxation curves of the experiments [|TTjl on glassy material with the remark- 
able accuracy show the algebraic decay rather than the standard exponential 
relaxation. The power-law relaxation can be expected to be a common feature 
of dynamical systems in a transition region between the stochastic and regular 
motion (supercooled liquids, glasses and polymer materials) jjDfl . 



3 Qualitative kinetic analysis of random processes with 
memory 



It is well known [|T3j that the generalized Langevin equation is of the form 



^ = - 1 J*M(t-t')v(t')dt' + L(t), 

where v is the particle's velocity, M some memory function, L the noise term. 
If the value v(t) is observed with a time resolution 5t ^> r c , where r c is the 
correlation time of forces producing the random particle motion, the Markov 
approximation is applicable, so that we obtain the ordinary Langevin equation 

In the case of complete but not ideal memory the generalized Langevin equation 
becomes a fractional differential equation. The evolution of the probability 
distribution associated with the velocity in the fractional Langevin equation is 
described by the fractional partial differential equation in the phase space for 
the phenomenon [[nj . 

Consider the probability distribution P(x, t) to find a particle at point x and 
time t. The normalization condition is 

DC 

J P(x,t)dx= 1. (10) 



- DC 



A broad class of various stochastic processes is given by the Chapman- 
Kolmogorov equation: 

= / M(x, t- x\ t') P(x f , t') #(f ) , (11) 

where the memory function M(x, t; x' : t') accounts for the probability distribu- 
tion P(x',t') in the previous instantes of time t' < t [lT5(| . The expansion of the 
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kernel M in terms of the difference x — x 1 yields the generalized (in the terms 
of the memory effects) Kramers-Moyal equation 

E("Vr/_* D^(x,t-t f )P(x,t')dt' : (12) 



dt 

where V = d/dx and the coefficients D^ 1 ' are the moments of the memory 
function M divided by n\. In the case of Markovian processes the moments are 
proportional to 5(t — t') and the integration in QT21) vanishes. It is relevant to 
remark that t in ([12]) should be treated in the "kinetic" but not "microscopic" 
sense. The Markov approximation (the Van Hove method |8j]) ignores memory 
effects in some sense, but such approach is not always useful. For example, the 
statistical theory describing transport properties of turbulent plasma leads to 
the conclusion that the turbulence is of subdiffusive nature and that the diffu- 
sivity considerably decreases. Therefore, memory effects can be important for 
explaining the dependence of the transport properties of saturated turbulence 
on the eigenfrequency of the unstable mode in the case of instability driven by 
the gradients in the coordinate space [jT^ . 
Using (|T2j) one can write the equation: 

r)P(r t) 00 1 r) n rt 

-^f 1 = E (-1)" -i j-,[Kn(x) I D(t - t>) P(x, t>) dt'] , (13) 

where K n (x) are arbitrary functions. The application of the Laplace transform 
to ([T3|) (with the initial values given on the whole real axis in the form P(x, 0)) 
leads to the following nonhomogeneous differential equation 

oo 1 3 n 

sP(x, s) - P(x, 0) = D(s) E ("If i jr- n [K n {x) P(x, s)] . 

n= i nl ox n 

Consider the case when memory is complete but not ideal (see previous para- 
graph). Then we come to the following equation 

sP(x, s) - P(x, 0) = § £ (-1)^ - ^~ n [K n {x) P(x, s)} , (14) 

s v n= i n\ ox n 

where D\ is the positive constant. Using the Caputo's definition ( |32| ) of the 
fractional derivative (see [[TTJ and Appendix) and its Laplace transform ([M|), 
the study of equation C P~5| ) for the case of complete but not ideal memory leads 
to the fractional generalization of the Kramers-Moyal equation: 

d 2(3 P(r t) 00 1 d n 
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If K\(x), K2(x) exist (nonzero) and K n {x) = for n > 3 the equation ( p75| ) is a 
fractional generalization of the Fokker-Planck-Kolmogorov equation | 



Observe some typical features of the fractional Fokker-Planck-Kolmogorov 
(FFPK) equation. If (3 = 1/2 [y = 0) it is transformed into the convential 
diffusion equation. This form is equivalent to the complete absence of memory. 
For (3 = 1 [y = 1) we have the convential wave equation, i. e. the process with 
complete memory. The equations containing derivatives of higher than second 
order with respect to time cannot exist in nature: a random process cannot 
spread faster than a collection of deterministic trajectories. (3 = defines the 
case of localization, which is the lowest limit of any diffusion process. Hence 
the physical bounds on (3 are given by < (3 < 1 (from ([T4|) it follows that 



1/2 < {(3= (z/+l)/2} < 1 because < v < 1). As is well known from statistical 
physics, one of simple criteria of irreversibility is whether or not equations are 
invariant in respect to time reversal (t — » —t). The specific character of the 
processes described by fractional time derivatives is that for the substitution 

{-tff 3 = ^{cos(2/?tt) + zsin(2/?7r)} 

the relative amount of system states is conserved, and the other one corresponds 
to irreversible losses |9j]. This allows us to suppose that for < (3 < 1 the FFPK 
equation describes the random processes with memory. 



4 Criterion for the relative degree of order 

The FFPK equation is a integro-differential equation in partial derivatives with 
varying coefficients, so in general one cannot find analytically its solution. For 
the stationary systems, the probability distribution does not depend on time 
(-^P(x : t) = 0), and their analysis leads to the Gibbs distribution, the corner- 
stone of statistical physics. A natural extension of the stationary analysis is the 
study of nonstationary systems in the self-similar regime, when the dependence 
on two arguments x, t is expressed in terms of a single variable y = x/a(t): 

P(x,t) = [a(t)] a ip(y), (16) 

where the functions a(t),(p(y) and the exponent a are to be defined. Mathe- 
matically, the probability distribution ([TBI) is a homogeneous function of order 
a. Physically, the transition to the new variable y = x/a corresponds to scaling 
the stochastic quantity x on an arbitrary scale a(t). It is well known[[n| that 



such a feature is displayed by fractal objects. If the domain of definition of 
the phase space for the stochastic system is a fractal set, its dimension D lies 
between 2 (the conventional phase space) and (the point of equilibrium). 
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In order to find the exponent a, we substitute the function ( |T6| ) into the 
normalization condition (|TU|), getting result 



a(t)}-^= J ip{y)dy. 



The left-hand side of this equation depends on the time, whereas the right-hand 
side does not. Hence, it follows that a = — 1. The form of the function tp{y) 
can only be found from the FFPK equation. 

Of different macroscopic functions, only the entropy S possesses a combina- 
tion of properties that allow to use it as a measure of uncertainty (chaos) in the 
statistical description of the processes in macroscopic systems [ 



Sit) = - J P(x, t) log(P(x, t)) dx + S = 

—oo 

oo 

= - / ip(y) log(p(j/)) dy + log(o(t)) + S . (17) 

— oo 

Define the mean value of x as 

oo 

(x 2 )=a 2 (t) J y 2 ip(y)dy = B-a 2 (t). 

— oo 

Using the above relation we can write the expression (|T7|) in the form 

oo 

S(t) = - J <p(y) log(tp(y)) dy - 0.5 log B + 0. 5 log(x 2 ) + S . (18) 



-oo 



The entropy S(t) depends on time only on account of (x 2 ). To consider the 
evolution of stationary states by means of slowly varying controlling parameters 
(they can be found among the parameters characterizing the stationary state) 
one can use the S-theorem [|20j as a criterion of the relative degree of order in 
various states which reveals for what states the degree of order is higher. It 
should be noted that the S-theorem considers only the stationary states using 
thermal equilibrium as the reference point for the degree of chaos. At t — > oo the 
probability distribution becomes negligible and it is necessary to renormalize 
the entropy: 

S = S{t)-0.5\og{x 2 ) . (19) 

The procedure is equivalent to fixing the value (x 2 ) for any value of a chosen 
controlling parameter (by the way, one of conditions of the S-theorem is the 
equality condition for the average effective Hamiltonian functions in different 
states of interest), since this influences only the reference point of entropy and 
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does not affect the renormalized entropy differences between different states. 
Although the value (x 2 ) depends on the controlling parameter (as well as on 
t), the fixing of the former does not mean fixing the latter. Now using the 
renormalized entropy difference as a measure of the relative degree of order, one 
can study the evolution of system states in the space of controlling parameters. 



5 Some Examples 

One of the simple cases of the FFPK equation is the Time Fractional Diffusion- 
Wave (TFDW) equation: 



Q2f3 Q2 

0^ P ^ = D 0^2 P ^t), (20) 

where D is a positive constant. Its fundamental solutions in case of the basic 
Cauchy and Signalling problems are well known [|2T|1. Let us take the point 



P(£, 0) = 5(£) as an initial position of a particle. If 1/2 < (3 < 1, it is necessary 
to specify the initial value of the first order time derivative QP(x, t)\t=o+, since 
in this case two linearly independent solutions are to be determined. To ensure 
the continuous dependence of our solution on the parameter (5 in the transition 
from 2(5 = 1~ to 2(5 = 1 + , we assume jLP{x,t)\t=Q+ = 0. Then the equation 
has the following solution 




P(x,t) = :r y^M[-^ ] f3}, (2i; 

where 

M(,,,3) = - L e—-^, O<0<1, 

Ha denotes the Hankel path (a contour that begins at a = — oo — ia (a > 0), 
encircles the branch cut that lies along the negative real axis, and ends up at 
a = -oo + zb (b > 0)). 

MainardifgTj developed the function M(z; (5) as a series 

oo (_ 1 \n n 



^ n\T[-/3n + (\-/3)\ 
and showed that it is a particular case of the Wright function 

W(z; A, M ) = E — tftv V = w~ e°+ z ° — 



where A > — 1, ji > 0|njj]. It is non- negative for any < ft < 1 and satisfies the 



normalization condition J n °° M((; (5) d( = 1. Clearly, these properties are also 
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characteristic to the probability distribution mentioned above. For (3 = 1 the 
function is the Dirac ^-function. For (3 = 1/2 we have the Gaussian function 
MO; 1/2) = exp(-z 2 /4)/ v / 7T- For < (3 < 1/2 the function M(z;(3) (z > 0) 
decreases monotonically, while for l/2</5<lit first increases and then 
decreases, exhibiting its maximum value at a certain point z max . The expression 
( gTj) describes the particle evolution in the space-time and is the Green function 
of the TFDW equation. Averaging x 2 we obtain 



* 2 )= J x2p ^t)dt= (22; 



r(2/? + i 

The probability distribution ( |2T|) obeys the following scaling relations 

P{x = b^x, i = bt) = b~P P(x, t) 

for an arbitrary parameter b. 

Note that if in (|2Tj) (3 = 1/2 we obtain the convential Brownian motion 
(purely random process). For (3=1 the particle performs ballistic motion 
(purely deterministic process). Next we are going to show that < (3 < 1 is 
the special case when chaos and ordered motion coexist. 

Now let us calculate the entropy (|T^) in this case: 

S(t) = 0.5 + 0.5 log 2 ((x 2 )) + S{(3) . (23) 

As a controlling parameter we choose (3 and will consider the evolution of the 
sequence of states corresponding to different values of the controlling parameter. 
One should keep in mind that lim^ 00 (^ 2 ) = oo. For this reason we must 
renormalize the expression ( P3f) according to flTDp. Then the relative degree of 
order can be estimated as 



3C 



S((3) = 0.5 log 2 T(2(3 + 1) - / M(z, (3) log 2 (M(z, (3)) dz. (24) 

o 

The dependence S((3) is represented on Fig. [[]. It reaches the maximum at 
(3 C = 1/2. We regard this state as a state of physical chaos (the correctness 
of this assumption will have to be verified). The smooth decrease of the value 
( ggp is a quantitative measure of the increase in the degree of order. Since 
the inequality 5(1/2) > S{(3 ^ 1/2) is satisfied, (3 —> h + A/3 is the transition 
from a less ordered state (physical chaos) to a more ordered state. This is an 
indication that we have found the corresponding controlling parameter, and the 
evolution of the system in the space of the controlling parameter is associated 
with self-organization. The conclusion is valid because for (3 = 1 we have 
the purely deterministic state which may be taken by a reference point of the 
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Figure 1: The dependence of the value S(/3) on the controlling parameter (3. 



degree of order. Some orderedness also happens to be the case for < (3 < 1/2. 
Although here we do not come to complete order, in some sense the orderedness 
is higher than for (3 C = 1/2. Thus the value (3 can be adopted as a measure 
of relative degree of order. It is useful to mention some particular values, 
lim^o S((3) = 1/ ln2 and £(1/2) = 0.5(1 + ln7r)/ ln2. To sum up, the process 
considered above is the special case when chaos and ordered motion coexist, 
the value f3 characterizing the relative degree of order in the process. Next we 
will compare this process with the fractional Brownian motion. 

In many physical systems the Gaussian distribution is a straightforward con- 
sequence of the Central Limit Theorem which makes it possible to consider com- 
pletely random processes. Hurst found a new set of statistical tools to examine 
the data which does not represent a purely random structure though standard 
statistical methods do not show any correlation between the observations^ . 
Using his rescaled range analysis one can extract meaningful information about 
the "memory" of a time-series. If the observations are not independent, each 
carries a "memory" of events which preceded it. Mandelbrot and Van Ness[j23| 
suggested the so-called fractional Brownian motion (the fractional integral of 
the white Gaussian noise) as a model reflecting the phenomenon. The fractional 
Brownian motion has the following probability distribution 

,2 \ 



Tl(x;t) = 



1 



x 



Its average (x 2 ) is proportional to t 2f} as in the case of the process (PT1). Then 



(25) 
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the expression (jig ) is of the form 

S(t) = 0.5{ln(27rD(f)) + 1} , (26) 

where D(t) = 2Kt 2/3 . Renormalizing it in accordance with ([191) we obtain 
S = 0.5(ln27r + 1). In other words, the renormalized entropy does not depend 
on (3, i. e. the fractional Brownian motions with various Hurst exponents are 
almost the same in respect to the relative degree of chaos. This example differs 
radically from the foregoing one: in the model of fractional Brownian motion 
( P5p the rise of orderedness with changing the Hurst exponent is impossible for 
the simple reason that the shape of the probability distribution does not change 
in the space of the controlling parameter (compare with the process ( PH) ). The 
numerical modelling verifies this conclusion j23j. 



6 Concluding remarks 

The fractional calculus formalism generalizing differentiation and integration to 
fractional orders has a long history [j215j], but recently the interpretative approach 



opened broad perspectives in physical and engineering applications |T2], 
The above consideration shows that fractional calculus provides a bridge be- 
tween purely deterministic processes and purely stochastic ones. The fact is of 
interest in its own right because chaos and order in Nature coexist. According 
to L. Boltzmann and J.W. Gibbs, in closed statistical systems, evolution in time 
results in the equilibrium state which is the most chaotic (or, in another words, 
purely random). The element motions in closed systems are independent of each 
other. In open systems the environment induces memory effects, so that the 
macroscopic behaviour of such systems contains a manifestation of microscopic 
dynamics. If microscopic motion in the systems is very complex or random, the 
complete memory transmits the complex behaviour to macroscopic evolution 
of the system as a whole. It is worth noticing that the model (used by us) of 
memory at the points of a Cantor set is not too exotic. It implies that memory 
is intrinsic to all time scales (in such a way that the corresponding memory 
function has no characteristic scales) of the phase space of a system given that 
the number of divisions generating a Cantor set tends to infinity. In a different 
way, the largest of time scales would be in the system as in the case of the 
Markov approximation. We have demonstrated that the relationship between 
the Cantor set and the fractional integral reduces the generalized (in the terms 
of the memory effects) Langevin, Kramers- Moyal, Fokker-Planck-Kolmogorov 
equations to their fractional form. The fractional generalizations turn out to be 
useful for studying the random processes with residual memory (without any 
characteristic time scale). 



6 CONCLUDING REMARKS 



14 



The process described by the Time Fractional Diffusion- Wave equation is 
an example clearly showing that for chaos and ordered motion in a system to 
coexist, its probability distribution must undergo qualitative changes in the 
function form with slowly varying the system parameters. In our consideration 
the parameter was the relative amount of the system states having ordered- 
ness during the system evolution. The evolution in the space of the parameter 
looks like a stochastic analog of bifurcation connected with the phase transition 
"order-disorder". For < (3 < 1/2 the probability distribution has one sharp 
maximum at z max = 0. If the parameter attains the value [3=1/2 the function 
becomes flat. Then for 1/2 < (3 < 1 the probability function takes the form 
with two maxima symmetrical in respect to the origin of coordinates (where the 
function minimum is found). When the parameter (3 goes towards 1, the prob- 
ability to find a particle between the maxima is becoming less and lesser, and 
the function peaks get narrower and higher. For (3=1 the probability distribu- 
tion is transformed in two (5-functions. As a result, this system state becomes 
completely ordered. The fractional Brownian motion has not the possibility. 
The Hurst exponent (which is the only parameter of the process) influences the 
asymtotic behaviour of its autocorrelation function, but the process remains 
Gaussian anyhow. 
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APPENDIX 

In this appendix we consider briefly the basic formulas used in the fractional 
calculus introduced in the text. The starting point is the introduction of the 
causal function $\{t) defined as 

$ A (£) : = | kxy *- 0; AeC, 

whose Laplace transform is 

— r+oo , I 

£[<S>x{t)] = $\{t) := / e- st ${t) dt = — , Re A > 0, Re s > 0. (27) 

jo § 

The function satisfies the composition rule 

<S>x(t)*<S>^(t)= f t <P x (r)<P^t-r)dr = <P x+ ^t), ReA > 0, Re/x > 0. 
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The integral of order n of a causal function f(t) can be expressed by the 
convolution between $ n and / 

basing on the well known formula that reduces the calculation of the n-fold 
primitive to a single integral. For A > 0, the function $>\(t) is locally absolutely 
integrable in < t < +oo. To extend the above formula from positive integer 
values of the index to any positive real values, let us define the fractional integral 
of order a > 0: 

I a f(t) := -Lj f Q {t - t)- 1 /(r) dr = * a (t) * f(t). (28) 

The Laplace transform of the fractional integral is the straigthforward general- 
ization of the ordinary case a = n 

C[I a f(t)] = % (29) 

For A < 0, the causal function $>\(t) is not locally absolutely integrable and 
consequently the integral 

Y^j£ jfi^ dT = *-(*)*/(*). (30) 

is in general devergent. Nevertheless, if A = — n, (n = 0, 1, . . .), the functions 
<&\(t) can treat in the framework of generalized functions In this case it 
reduces to the n-derivative (in the generalized sence) of Dirac (5-function 

Then formally the derivative of order n of a causal function f(t) can be obtained 
by the convolution between 4>_ n (£) and / 

—fit) := J o f(r) 5^(t -r)dr = *_ n (t) * /(*), t > 0. 

The limit case a = defines the Identity operator 

ft 

In order to obtain a definition for the fractional derivative (with non-integer 
a) that is valid for classical functions, we have to regularize the divergent in- 
tegral (BO) in some way. As a consequence we arrive at two alternative regular 



]*f(r)6(t-T)dT = $o(t)*f(t) = f(t). 
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definitions for the fractional derivative D a , which read for n — 1 < a < n 

D cf(t) = IT? — ~ f , f{ \\, dr, (31) 

c J w r(n - a) dt n Jo (t - T ) a+l - n v ' 

d a 1 rt f^(r) 

Drf(t) = = tv~ ~T /„ ' L,-, dr. (32) 



eft" v 7 r(n - a) Jo (t - r) a+l 
The difference between the two definitions lies in 

71—1 

D a J{t) = D°f(t) + E / (fc) (0 + ) $(*-«+!)(*). 

The defintion (^T|) is the most commonly adopted in mathematically oriented 
papers and books []261|. The definition (|520 , introduced by Caputo [|T7fl, is more 



restrictive than flgTp because it requires the function /(£) to be n-defferentiable. 
As it was shown in this paper, the latter definitions is more suitable for the 
problems considered above where the conventional initial conditions are ex- 
pressed in terms of integer derivatives. Using the classical technique of Laplace 
transform 

An n-l 

C[ h f{t)]= S " /(S) ~ E sn " 1 "V ( * ) (0 + ), n e N (33) 
and (|27]), fl32j) we get 



Aol n -i 

£l^f(t)}=s a f(s)-Y,s a ~ 1 - k f {k X0 + ), n-l<a<n, (34) 
at n= Q 



first stated by Caputo [|T7j. It is worth noticing that according the definition 
( ftll) the fractional derivative of a constant does not vanish if a is not integer, 
while according to the defintion (1321) it vanish for any a. 
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